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We study the character of an Ising nematic quantum phase transition (QPT) deep inside a d- 
wave superconducting state with nodal quasiparticles in a two-dimensional tetragonal crystal. We 
find that, within a 1/N expansion, the transition is continuous. To leading order in 1/A'^, quantum 
fluctuations enhance the dispersion anisotropy of the nodal excitations, and cause strong scattering 
which critically broadens the quasiparticle (qp) peaks in the spectral function, except in a narrow 
wedge in momentum space near the Fermi surface where the qp's remain sharp. We also consider 
the possible existence of a nematic glass phase in the presence of weak disorder. Some possible 
implications for cuprate physics are also discussed. 



I. INTRODUCTION 

In this paper we study the nematic to isotropic quan- 
tum phase transition (QPT), deep within the d-wave su- 
perconducting phase of a quasi two dimensional tetrago- 
nal crystal. The "nematic" refers to a broken symmetry 
phase in which the four-fold rotational symmetry of the 
crystal is broken down to a two- fold symmetry. More 
specifically, it is an "Ising nematic," which spontaneously 
breaks the discrete rotational symmetry of the tetragonal 
crystal to an orthorhombic subgroup, — > C2v, while 
retaining translational symmetry. The superconducting 
order opens a gap in the quasiparticle excitation spec- 
trum, except at four gapless nodal points, but these nodal 
quasiparticles (qp's) couple strongly to the nematic or- 
der parameter fluctuations. In the nematic phase, these 
nodes are displaced from the relevant symmetry direc- 
tions by an amount proportional to the nematic order 
parameter, as shown in Fig. [TJWe derive a phenomeno- 
logical theory to describe this "nodal nematic QPT" con- 
sisting of a nematic mode coupled to the nodal qp's of 
the d-wave superconductor. 

Our motivation to investigate this problem is two-fold; 
Firstly, there is, now, considerable experimental evidence 
that a nodal nematic phase occurs in at least some "un- 
derdoped" cuprate superconductors, so this study has 
potential relevance to the transition from this state to 
the isotropic state in these materials. The best ev- 
idence of this comes from measurementJ^ of strongly 
temperature dependent transport anisotropics in un- 
derdoped YBa2Cu306-H5, and more recent (and more 
direct) neutron scattering experiments in underdoped 
YBa2Cu306.45 (YBC0)2l. Specifically, the spontaneous 
onset of one-dimensional incommensurate spin modula- 
tions in the neutron experiments is clear evidence of an 
isotropic-to-nematic transition with transition tempera- 
ture Tn ~ IhQK. For 5 = 0.45, is greater than 
the superconducting Tc — 2>bK. However at higher O 
concentration, 5 > 0.7, both the neutron scattering and 
transport experiments see no evidence of a transition to 



a nematic phase down to the lowest temperature d^ * '^ ! "^ -*^. 
It is thus reasonable to assume that Tn{S) ^ for a 
critical value, 0.45 < Sc < 0. 7, with a QPT at 5 = (5c 
inside the SC phase. l^^lSMlSHI]. The extent to which ne- 
matic phases are generic to the cuprates is a question 
that is beyond the scope of the present study. However, 
Matsuda and coworkers^^ have recently reported that the 
"fluctuating stripe" phase of underdoped La2_a;Sra;Cu04 
(with X — 0.04) is act ually a (diagonal) nematic phase. 
Moreover, STM studied^^^ have revealed a glassy phase 
with nematic domains in underdoped Bi2Sr2CaCu208+5. 

Secondly, the topic of quantum critical points (QCP) 
in 2D systems with itinerant fermions has been and con- 
tinues to be a topic of broad interest. Here we study 
a new QCP in a 2D itinerant fermion system. In gen- 
eral, the presence of gapless fermions at the Fermi sur- 
face makes the study of this topic theoretically challeng- 
ing. Since gapless fermions interact with massless bosons 
associated with order parameter fluctuations, this inter- 




FIG. 1; Quantum critical point at A = Ac, separating the 
nodal nematic phase for A^^ < A^T^ from the symmetric 
phase, and its quantum critical fan. T^ and r„ are the super- 
conducting and nematic critical temperatures and the (pur- 
ple) wedge corresponds to the thermal critical regime. 



action affects both low energy degrees of freedom, making 
it difficult to obtain a correct description of the critical 
physics. Simplifications occur when one considers nodal 
fermions which are gapless only at four nodal points of 
the underlying Fermi surface in a d-wave superconductor. 
The limited phase space for the gapless Fermions restricts 
the possibilities for scattering mechanisms permitting a 
controlled analysis. 

Electronic nematic phases were first predicted to occur 
in doped Mott insulators"^, and have now been experi- 
mentally observed in a number of different systems. In 
addition to YBCO, they have been discovered in trans- 
port experiments in semiconductor heterostructured^ 
and in the bulk transition metal oxide Sr3Ru207ESl. 
While in all cases the thermal transition to the nematic 
phase appears to be continuous, the character of the 
quantum transition at zero temperature is less clear, and 
at least in the case of Sr3Ru207'^ appears to be first 
order. However, interesting universal physics can be ex- 
pected in the case of a continuous transition. 

Vojta, Zhang and Sachdev'^'^^' analyzed various types 
of ordering transitions in a nodal superconductor, includ- 
ing Ising nematic ordering. They employed an e expan- 
sion (e — 3 — d) to derive a perturbative renormalization 
group (RG) flows near the decoupled fixed point, A = 0. 
In the nematic case, they found runaway flows which they 
tentatively interpreted as a fluctuation induced first or- 
der transition. Of course, a runaway flow can also imply 
a breakdown of perturbation theory. 

In this paper we work directly in d = 2, and use a 
large N theory to access a non-trivial critical point at 
finite coupling Ac. This is achieved by generalizing the 
problem to one in which there are A'' "flavors" of nodal 
qp's (with the physical N = 2 case corresponding to two 
spin polarizations/^. In the large N limit, we present a 
theory of the nodal nematic QPT^^, whose critical behav- 
ior should be smo othly and systematically corrected by 
the 1/iV expansioEpMmil], We find: a) the quantum 
phase transition is continuous, with non-trivial critical 
exponents that we compute for iV — s- oo, b) the transi- 
tion occurs at a finite critical coupling between the nodal 
qp's and the nematic order parameter, so is inaccessible 
by perturbation theory, c) the anisotropy of the nodal qp 
dispersion strongly influences the interplay between the 
critical fluctuations and the qp's. We also speculate on 
the implications of this theory to the phenomenology of 
nodal qp's in cuprates. With minor changes, this the- 
ory also applies to a possible electronic hexatic QPT in 
graphene and similar systems 

The paper is organized as follows. In Section [IT] we 
introduce the model and discuss its symmetries. The ef- 
fective theory of the nematic order parameter is derived, 
in the large N limit, in Section |III[ In Section |IV| we 
discuss the behavior of the nodal quasiparticles at and 
near the nematic quantum critical point. In Section [V] 
we give a qualitative presentation of of the physics of a 
nematic glass and how it affects the spectrum of nodal 
quasiparticle, leading to a Fermi arc phenomenology in 



the quasiparticle spectral function. In Section VI we dis- 
cuss a number of important open problems. 



II. THE MODEL 

The effective Lagrangian which describes the coupling 
of the nodal fermionic qp's to the nematic order param- 
eter is C = + £int + J^4>, where is the linearized 
nodal qp Lagrangian in a pure dx^~y^ SC 



Here 



(2.1) 



(2.2) 



are two-component Nambu spinors representing the 
nodal qp's for node index n — 1,2 and a, f3 = I, . . . ,N are 
Sp{N/2) "flavor" indices which, for N = 2, correspond to 
the qp spin polarizations, since 5*^(1) ~ SU{2). Eq.(2.2l 



represents nodal qp's at momentum positions p relative 
to each pair of nodes at Ki = {K, K) and K2 = {—K, K) 
for n = 1 , 2 respectively and 
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(2.3) 



are velocities normal and tangential to the Fermi surface 
(FS). Note that we investigate our model for general val- 
ues oivp and i;a- We use n and T3 to denote 2x2 Pauli 
matrices acting on the Nambu spinors. The Lagrangian 
for the nematic order parameter (j) is 
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u 
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(2.4) 



where the ellipsis represent higher order terms in powers 
of the nematic order parameter field 4>. We assume that 
u > and that > 0, i.e. in the absence of coupling to 
the nodal qp's, the system is in its isotropic phase, {(/)) = 
0, and the nematic mode is gapped. The interaction term 
that couples the nematic order parameter to a fermion 
bilinear is 



(2.5) 



which has a well defined large N limiiPI^SIlIIII]. This 
model has the discrete symmetry (0 — > — ^ 
—p, ^' T2^'). The combination of this discrete sym- 
metry with the implicit symmetry of the Nambu nota- 
tion and the two fold rotation symmetry associated with 
{"^i ^ ^'2, kx ^ ky) amounts to symmetry. Should 
{(j>) ^ 0, the position of the nodes shift in the v^, direction 
and the 6*4^ symmetry is spontaneously broken down to 
C2V associated with (vPi ^ ^'2, kx ^ ky) (see Fig. [T]). 
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The interaction in Eq.(2.5l has the form of a Dirac The fermion determinant in Eq.(3.2l represents the sum 



fermion current coupled to one component of a vector 
gauge field. Couplings of this type have been studied 
in the context of gauge theoretic approaches to nodal 
fermions and phase fluctuations^^ However, in the 

case at hand, the microscopic theory does not have a local 
gauge invariance associated with the nematic field 0, but 
instead only has a global discrete symmetry which can be 
broken spontaneously. 

Since the coupling makes the tangential and normal 
directions to the FS inequivalent, one must consider an 
anisotropic dispersion, vf/va 7^ 1- Indeed, in previous 
studies of quantum criticality with Dirac-like fermions, 
an effective rotational symmetry about the nodal point, 
''^f/va ~* 1; was shown to be an emergent property of the 
fixed point in the sense that anisotropy is perturbatively 
irrelevant; this is not true in the present problem, and 
indeed we find a number of qualitatively new features of 
the critical theory which derive from this anisotropy! 

The fact that vp and wa have different physical ori- 
gins, implies that there is no reason a priori to expect 
the existence of even an approximate symmetry vp — va- 
Indeed, experimentally it is well known that the degree 
of anisotropy is extreme: vf/va 19 for BSCCO from 
thermal conductivitjI^H and ARPE^, and vf/va ^ 14 
for YBCO from thermal conductivity^. After we formu- 
late our theory for general values oi vp and va, we inves- 
tigate the specific case which corresponds to linearizing 
the phenomenological band structure of M. Norman et 

aim 



III. THE ORDER PARAMETER THEORY 

We will now proceed to study the behavior in the large 
A'' limit of the theory of nodal fermions coupled to a ne- 
matic order parameter theory presented in the last sec- 
tion. As it is standard in the large A^ limit we proceed 
to formally integrate out the fermions and obtain the re- 
sulting effective action for the nematic order parameter 
ficld.^^'^^^^ In the large-A^ limit, the resulting bosonic 
path integral 
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(3.1) 



will be dominated by the saddle point solution to the 
effective action S'off [0], and is amenable to a controlled 
1/A expansion to evaluate the effects of fluctuations. In 
imaginary time, the full effective action S'off[(^] has the 
standard form 
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of fermion bubble Feynman diagrams. Since the action 
of the nodal fermions has the structure of a anisotropic 
Dirac theory with a coupling to the nematic order pa- 
rameter (j), each bubble diagram involves a trace over a 
product of Dirac fermion propagators with the appropri- 
ate insertions of the matrices ti dictated by the form of 
the coupling. This is a standard procedure which ob- 
viates the need to present the well known intermediate 
steps. We will thus only present the main results of these 
calculations. 

Upon a rescaling the order parameter <j) ~^ 0, it is 
apparent that the effective action for the rescaled fleld 
has the form NScs[4']- Hence, in the large A'' limit the 
path integral is dominated by the solution to the saddle 
point equations for the effective action S'off[0], and its 
fluctuations. When expanded in powers of 4>, the contri- 
bution of the fermion determinant introduces a negative 
contribution to the quadratic ter m that is proportional 
to in the effective action [</>] .HaiMlil This allows for 
the spontaneous symmetry breaking dri ven by t he cou- 
pling constant A. In theories of this typ^i^IIlEl] ^j^gj-g jg 
a critical coupling strength Ac, beyond which 4> gains a 
finite expectation value due to the spontaneous symme- 
try breaking. The critical value of the coupling constant, 
the location of the quantum phase transition, depends on 
the bare mass of the nematic mode and the UV cutoff. 
If we use a sharp UV cutoff A, in the large N limit we 
obtain 



Xl = [2TT^/ln2]{m'^VFVA/A). 



(3.3) 



For A^ < Ac the symmetry is unbroken and {(p) = 0. 
However for A^ > A^, a nonzero expectation value {(j>) is 
obtained by minimizing ^cff This can be achieved by 
solving the saddle point equation: 



'■{<!>) = - 




For A close to Ac, the solution of Eq.(3.4l satisfies 
(0)oc(A-Ac)i/^ 



(3.4) 



(3.5) 
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Thus, in the large A^ limit, we find a continuous quantum 
phase transition, i.e. a quantum critical point, separat- 
ing the weak coupling isotropic phase at A < Ac from a 
strong coupling phase at A > Ac. Th e critical exponent 
of the order parameter of Eq. (3.5), 1/2, obtained in 
the large A^ limit will be corrected order by order in the 
1/A expansion. Within the 1/A^ expansion, the cor- 
rections to the value of this and other critical exponents 
of this theory follow from the resummation (exponentia- 
tion) of the infrared divergent 1/A^ perturbation theory. 
This procedure can also be viewed as a construction of 
a fixed-point theory (and its associated renormalization 
group) order-by-order in the 1/A^ expansion, a standard 
procedur e in relativistically invariant theories of Dirac 
.2) fermions.lialll] 
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In the present problem, Lorentz invariance is explicitly 
absent due both to the velocity anisotropy, parametrized 
by an anisotropy ratio a = [vp — v/s^/vp, and the form of 
the coupling between the nematic order parameter </> and 
the Dirac fermions. The lack of this invariance changes 
the physics in dramatic ways. As noted above, Vojta, 
Zhang and Sachdev studied the quantum critical behav- 
ior of a similar system by means of a perturbative renor- 
malization group (RG) analysis at fixed N near the crit- 
ical (renormalizable) dimension Z? = 3 -I- 1 using an e 
expansion.^ Due to the combined effects of the veloc- 
ity anisotropy and the anisotropic coupling, they found 
that the perturbative RG has runaway flows, signaling 
the breakdown of perturbation theory. Following a stan- 
dard analysis (see for instance Refs. |21l23j ) these au- 
thors interpreted this runaway flow as an indication of a 
fluctuation induced first order transition. Instead, in the 
large- iV theory that we present here, we find that, at fixed 
dimension D = 2 + 1 but for very large N , the quantum 
phase transition is continuous and has a finite renormal- 
ized velocity anisotropy which plays a key role.l^s! This 
behavior also contrasts with that of more nearly Lorenz 
invariant models, such as the QED3 type models studied 
in Refs. |25I26I27] . where a small velocity anisotropy was 
found to be either irrelevant or a redundant. 

We will now proceed with the analysis of this theory in 
the large- iV limit. The leading quantum fiuctuations at 
large N are obtained by expanding the effective action to 
quadratic order about the saddle point (see for instance, 
Ref. [31]). This standard procedure, when apphed to 
the present problem which breaks Lorentz invariance, in- 
troduces two separate energy scales associated with any 
given momentum p 



term generated by integrating out the nodal fermions. 
The scaling dimension of the nematic field (f, as tradi- 
tionally defined, is dim[(/3] = (d + z — 2 + 77)72, with 
the space dimension d, the dynamic critical exponent z, 
and the anomalous dimension 77. In the present case, 
dim[(^] = 1 and z = 1 is inherited from the Dirac- like 
spectrum of the fermions. This can then be interpreted 
as a (large!) anomalous dimension rj = ll^^l. As a result 
of this scaling, all local interaction terms of higher order 
in ip are (sometimes dangerously) irrelevant operators. 
(This is in contrast to a local ip'^ theory in which the wp'^ 
term is relevant for d < 3 at the Gaussian fixed point, so 
long as z < 2, making the fixed point unstable.) Non- 
local interactions are also generated by integrating out 
the fermions, which, presumably, lead to 0{\/N) correc- 
tions to the various critical exponents. Away from criti- 
cality, the correlation length scales as ^ ~ |Ac — Aj^'^ with 
V ~ \. To the extent that scaling holds, Tc ^ {\ — Ac) 
for A > Ac {vz — 1). The resulting phase diagram, with 
a V-shaped quantum critical fan is sketched in Fig. [T] 
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where Ei^px^Py) — E2{py,Px)- In terms of Ei{p) and 
E2 (p) , the Gaussian action S*-^' [tp] for the fluctuations of 
the nematic mode p{x) deflned by ip{x) = j^[(j){x) — {(j))] 
takes the following form on the disordered side (A < Ac) 
of the critical point: 



(Ppdu! 1 



(27r)3 2 
i^uj^ + Ei{p) (1 



{Xl~X')+Kico' + c^p') 

[Px ^ Py] 



viPy 



LU^+E^ip) 



(3.7) 



where k = Xcjrr? and 7 = y^^jZlvpVA- The quartic and 
higher order terms in fluctuations have coefficients which 
are^^ 22 (Jq-^^^ powers in 1 jN . Here we will only treat 
the physics to leading order in this expansiorP2l 

We can analyze the effective action in Eq. (3.7 1 from 
the viewpoint of scaling theory. Clearly, at criticality 
(A = Ac), the bare dynamical terms in Eq.(2.4l are irrel- 
evant, as they have larger dimension than the non-local 
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FIG. 2: The nematic mode spectral function B{ui,p) plot- 
ted as a function of energy lj/A at momentum p = 
(0.25A/uf,0.05A/uf), where A is a UV cutoff and is the 
Fermi velocity, (a) Isotropic case vf ~va~v. (b) Anisotropic 
case vf/va = 19.5. 



We compute the spectral response of the nematic 
mode by an analytic continuation of the effective action 
Eq. (3.7) to real time. The irrelevant terms in the effec- 



tive action can be set to zero for energies small compared 
to 7/k. In Fig. [2] we show the spectral function of the 
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nematic mode at criticality 

B{p,uj) = -2sgn(w)ImG'(p,w), 
for the if propagator 



(3.8) 



G(p,a;) = [^^-oj^ + E^^VM^) + (1 <"> 2)] \ 

(3.9) 

where we defined projectors 
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The spectral function of Eq.(3 



(3.10) 

contains contribu- 
tions from a pole and from branch cuts. When the disper- 
sion is isotropic, the nematic mode pole is at w = v\p\/^/2 
with residue Z — 1/4 and there is a threshold at w = dp| 
for the continuum due to the branch cut. (See Figj2][a/^ 
However, the dispersion anisotropy a > {vp > va) 
completely alters the analytic structure of the nematic 
mode spectral function. First of all, the pole approaches 
the continuum with a reduced residue^'*. Moreover the 
threshold energy scale for the continuum ai uj — v\p\ 
associated with given momentum p for an isotropic dis- 
persion, splits into features at two energy scales Ei{p) 
and E2{p) (See Fig|2jb)). The spectral function B{q,uj) 
can in principle be measured with momentum resolved 
inelastic X-ray scattering. Observation of nematic mode 
spectral function with two distinct energy scales associ- 
ated with a given momentum could serve as direct evi- 
dence for nodal nematic criticality. Moreover, the non- 
trivial analytic structure of the nematic mode spectral 
function results in momentum dependent scattering for 
quasi particles. 



IV. PROPERTIES OF NODAL FERMIONS AT 
CRITICALITY 

The critical nematic fluctuations have drastic effects 
on the nature of nodal fermions. To illustrate this, here 
we will discuss the single particle spectral function at the 
nodal nematic QCP, defined in the vicinity of each node 
n = 1, 2 as 



An{p,u) = -2sgn(w)Im[t/„,ii(p, w)]. 



(4.1) 



Here Gn,ii{p,^) is the (11) component of the (2 x 2) 
Nambu matrix time ordered qp propagator 



gnip,uj) = -i(r^'„(-p,w)«'„(p,w)). 



(4.2) 



Note, in the rest of this section, we will focus on the 
vicinity of the n = 1 node and omit the subscript n = I 
to simplify the notation. 

A matrix- valued qp self-energy S(p, w) best charac- 
terizes the effect of critical nematic fluctuations on the 
single fermion spectral function. We can then express 
A{p, to) in terms of S in the standard manner. In the 



presence of S, the Nambu matrix propagator for the 
nodal Fermions Q is given by Q^^ = Qq^ — S, where 
Qq{p,u!) ~ (luI — vppxT^ — VAPyTi)~^ IS the Nambu ma- 
trix propagator of the free nodal Fermion theory. Decom- 
posing the matrix valued self-energy in a basis of Pauli 
matrices as S(p,cj) = - S^^Vg + S^^Vi, the asso- 

ciated single particle spectral function is 



A{p, uj) 



Im 



~2sgn(aj)x 

(c.-S(o)) + («;^p,-SW) 



(4.3) 



From Eq.(4.3l one can understand the effect of differ- 



ent components of S on the physical properties of the 
fermions. 

To order 1/iV, the (2 x 2) qp self-energy matrix at the 
QCP is 



S(p,^)=; 



2N J (27r)2 27r 



Ti^o(fc, oj')TiG{p+k, u+uj') 
(4.4) 

where G{k,u}) is the large N nematic mode propagator 
of Eq. ( 3.9[ ). Notice that S'^*)(p', w) have explicit non- 
trivial dependences on p and uj. This momentum de- 
pendence makes it challenging to obtain the self-energy 
in closed form in general. Hence we will evaluate the 
self-energy and the resulting spectral function numeri- 
cally. Nonetheless some valuable analytic understanding 
can be obtained in limiting cases. First of all, in the limit 
of p = (nodal point), we can show that the self-energy 
takes the forrrP^ 



E(p = 0,cj) oc -^/(a) 



In 



A2 



(4.5) 



where f{a) is a function of dispersion anisotropy ratio a. 
Thus the imaginary part of self energy grows with energy 
and the width of the quasiparticle peak is never sharper 
than the energy. 

Another important observation we can make from 



Eqs.(4.4 - 4.51 is that the real part of self energy 
renormalizes va downwards without affecting vp- This 
is a combined result of two aspects of the theory. 
For one thing, each component I]*^*-'(p, w) scales with 
UJ, vfPx, VAPy respectively for i = 0, 1, 2 as it can be eas- 



ily checked from (4.5 1. Secondly, S'^^) enters the spectral 
function with the opposite sign from S^*^^ and S^^^ due 
to the structure of the interaction vertex, which reflects 
the fact that the coupling to the nematic mode breaks 
Lorentz invariance. Since there is a logarithmic diver- 
gence in this renormalization of va, the anisotropy ratio 
a is relevant at a = 0. Hence it is clear that one should 
consider the case of a anisotropic dispersion a =/= from 
the beginning. 

It is noteworthy that for non-zero p, any amount of 
anisotropy a ^ qualitatively changes the analytic struc- 
ture of the nodal Fermion self-energy S(p, w) from that 
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FIG. 3; Momentum distribution of the nodal nematic QCP spectral function at a energy ijj — —9meV for =0.508eV(7r/a) ^, 
VA = 0.026el/(7r/a)~^, {vf/v/\ = 19.5). Here, the momentum p is measured with respect to nodal point Ki = {K,K) in a 
(rotated) local coordinate system so that px and py lie, respectively, along the the nodal and the tangential directions, (a) 
Contour plot (color online). 8c is the critical angle at which a well defined qp peak departs from the incoherent continuum. 
For d < 6c, the qp is well defined, (b) and (c) are line cuts along the nodal direction and tangential direction respectively. 



at a = 0. This is due to the existence of two energy 
scales Ei{p) ^ E-iip) (see Eq.(3.6l for the definition of 
these energy scales) entering the nematic mode propa- 
gator Eq.(3.9l (see Fig. |2|. Such change in the analytic 
structure is not perturbatively accessible in a from the 
a = case. Hence we have computed S numerically 
from Eq. (4.4 1 for an arbitrary nonzero value of a, by 



obtaining ImS via a Monte Carlo integration, and from 
this obtained ReS by Kramers-Kronig. We verified this 
method against the analytic expression Eq.(4.5l in the 
isotropic limit. 

We present the nodal nematic QCP fermion spectral 
function A(j), cj), obtained from the numerical calculation 
of S(p, a;), in Fig [s] and Figji] Although we have car- 
ried out the calculation for arbitrary values of velocities, 
we used the dispersion vp — 0.508ey(7r/a)^^ and va — 
0.026eT^(7r/a)~^ with the velocity ratio vp/vA = 19.5 in 
our plots in order the demonstrate the effect of a large 
bare velocity anisotropy. This dispersion was obtained 
by linearizing the phenomenological model Hamiltonian 
for BSCCO of M. Norman et alPl. For the value of the 
critical coupling, we used \\/2N^ — 0.3. Fig [s] shows 
the momentum distribution of the spectral function at 
a fixed energy oj = —9meV. The eccentricity of the el- 



lipses in the contour plot refiects the large bare velocity 
anisotropy. We also show two line cuts or momentum dis- 
tribution curves (MDC's): one along the nodal direction 
and the other along the tangential direction. In Fig|4]we 
show representative energy distribution curves (EDC's) 
at two fixed points in momentum space (one point along 
the nodal direction, another point along the tangential 
direction) . 

On the basis of this numerical calculation of E(p, a;) 
and hence the spectral function A(p,u;), we extract two 
principal qualitative effects of the nematic critical fluctu- 
ations on the nodal Fermion properties: (i) strongly mo- 
mentum (angle) dependent scattering; (ii) an anisotropic 
renormalization of the already anisotropic bare velocities. 

The presence of strongly angle-dependent scattering is 
clearly shown in Figjsjja). The qp's are highly damped 
in the direction normal to the FS {px direction) while 
they remain sharply defined in the vicinity of the FS {py 
direction). More specifically, the momentum distribu- 
tion of the spectral intensity exhibits distinct behaviors 
in different wedges of p-space, which open about the FS 
at a critical angle, dc — ta,n~'^{vA/vF). For 6 > 9c, the 
peaks in the spectral function are broad with a width 
cx For 9 < 9c, the MDC has sharp peaks (the thin 
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FIG. 4: EDO's taken at two points in momentum space, (a) 
At p = (-9me\//uF,0). (b) At p = (0, 9me\//uA). These 
are the momenta corresponding to the position of the peak in 
Fig. 3(b) and to the threshold in Fig. 3(c) respectively. 



anisotropy. Such anisotropic renormalization of the dis- 
persion is the result of the structure of the coupling be- 
tween the nematic mode and the nodal fermions and it 
holds order by order in the (1/iV) expansion. The con- 
sequences of a RG analysis of this decrease in wa will be 
presented in a forthcoming papei^^^. 

One can readily understand the kinematic origins of 
the sharply defined qp's near the FS inside the narrow 
wedge in Fig. |3] (a). Ei{'p) defines the bare dispersion 
of a qp near the nodal point Ki = {K,K) while i?2(p) 
defines the bare dispersion of a qp near the nodal point 
K2 — {K,—K). Since the coupling to qp's near both 
nodes determines the dynamics of the nematic mode, it 
is the lesser of Ei (p) and E2 (p) that sets the threshold for 
decay. For p= {px,0) ~ {—9meV/vF,0) normal to the 
FS (Fig.ga)), E2{p) = va\p.\ <^ Ei{p) = vf\Px\, so the 
qp is highly damped. Notice the asymmetric line shape 
which reflects the p^dependence of the self-energy. For 
p= (O^py) = (0,9mel//uA) along the FS (Fig. gb)), the 
nematic fluctuations renormalize the qp velocity ua — 
VA, and the qp energy, Ei{p) = v^lPyl < £^i(p) < E2{p), 
so there no damping. It also produces an EDC with a 
'peak-dip-hump'. The existence of entirely undamped 
quasiparticles inside sharply defined fc-space wedges is 
likely an artifact of the first order corrections in the 1 /N 
expansion, since higher orders terms are likely to intro- 
duce finite damping. It is therefore more plausible that 
the wedge delineates a crossover from a regime in which 
the qp's are highly damped, to a regime inside the wedge 
in which the qp peaks are relatively narrow with a width 
that vanishes as one approaches the tangential direction 
along the FS. 



white lines) dispersing with a rcnormalized velocity va, 
enhancing the eccentricity of the constant energy ellipse. 
We have marked the boundary between different wedges 
with the dashed (pink) lines in Figjsj^a). Due to the ex- 
treme bare anisotropy, the angle 29c which defines the 
region in p-spa,ce with well defined qp's is rather small. 
However, this narrow wedge can qualitatively affect the 
long time properties of the nodal fermionJ^. 

The anisotropic renormalization of the dispersion is 
particularly evident in the Figjs] (b) and (c). The free 
nodal fermion theory would have placed sharply defined 
peaks at p— {\w\/vf,0) and p = (0, |w|/fA), in Fig|3](b) 
and (c) respectively. Figjsjb) indeed shows a peak (al- 
beit broad) at the position expected from the bare value 
of the vf and hence vf is unrenormalized by 1/iV fiuctu- 
ations. However, the position of a sharp peak in Fig|3jc) 
has been shifted away from the bare position. This shift 
is characterized by a renormalization of va ^ va, where 
■Da is related to va by 



N \vaJ 



O 



iV2 



(4.6) 



Here T{vf/va) is a positive function of vp/vA- Hence 
critical fiuctuations effectively enhance the dispersion 



V. THE NEMATIC GLASS 

Consider now the nodal nematic phase away from crit- 
icality, A > Ac, where the nematic order parameter has 
a non-zero expectation value, as in Eq.(3.5l. Because 



of the Ising character of the ordered state, there is a 
gap in the nodal nematic fiuctuation spectrum, wjv ~ 
T/v ~ (Ac — A). At energies large compared to wjv, the 
spectral function is little different from its behavior at 
criticality. However, at energies small compared to ujn, 
the nodal qp's are undamped, but the nodal positions 
are shifted from the symmetric points in fc-space by an 
amount Ap cx \/A — Ac as shown in Fig. ^ and the ve- 
locity tang ential to the Fermi surface is renormalized, as 
in Eq.(|46l). 



Unfortunately, quenched disorder (impurities) has a 
devastating effect on nematic phases, and indeed macro- 
scopic manifestations of electron nematic order have only 
been seen in ultra-pure systcms^^ Specifically, the 
order parameter theory in the presence of impu rities 
is equivalent to the random-field Ising modeP^^, and 
hence nematic systems are more generally expected to ex- 
hibit glassy dynamics rather than a broken symmetry.l^Sl 
Thus, more detailed spectroscopies and/or local mea- 
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surements play an especially central role in experimental 
studies of nematic phases. Here we sketch the effects of 
disorder on the single-particle spectral function in the 
nematic phases. We consider two extremal cases: short 
range and mesoscopic range disorder. 

Short range disorder produces an imaginary part to 
the self-energy, which smears the MDC at a; = in all 
directions about the nodal point. Due to the intrinsic 
dispersion anisotropy, this broadened weight is severely 
elongated along the tangential direction so the resulting 
MDC has a shape that resembles a Fermi "arc." However 
the width of the MDC in the normal direction is still 
proportional to the arc length. 

In sufficiently clean systems, the effects of disorder are 
only significant on longer distances, producing a "ne- 
matic glass" . Here we focus on the static effects leav- 
ing out (important) dynamical effects generally expected 
in glassy phases, such as hysteresis and noise. If each 
domain is large enough to support a defined local ne- 
matic order, the variance = among domains 
amounts to a distribution of nodal positions along the 
FS. In such a mesoscale glass phase, a local measurement 
(such as STM) will see sharply defined nodal quasiparti- 
cles at low enough energy, but a measurement, such as 
ARPES, which averages over a large area, will measure 
the domain averaged spectral function 

A{p,uj) ^ Jd(j)P{(t)) Aoip^, Py + (f>/vA, uj) (5.1) 

where -P(0) is the probability density that (f> which 
shifts the position of the node by Apy = (fi/vA, and 
AQ{px,Py,uj) is the spectral function for free qp's 



AoiPx,Py,uj) 



6{uj - Ei{p)){Ei{p) + VFP.) 



Ei{p) 

+ 5iLu + Eiip)){Ei{p) - VFPo.) (5.2) 



with Ei(p) given in Eq. (3.6 1. For a Gaussian distribu- 
tion of (j) with variance a, P{4>) — e"*^ Z^*^ /-v/Sttct and 
it is straight forward to show that the average spectral 
function in Eq.(5.1l becomes 



^Up' - {vFPxf 

e ^ -|- e 



(5.3) 

Note that this defines a Fermi arc of length (t/va with a 
vanishing perpendicular width, along the FS (j)^ = 0) at 
zero energy w = 0. 



VI. DISCUSSION 

In this paper we have presented a phenomenological 
theory of the nodal nematic QPT. Formally similar prob- 
lems arise in studies of the nodal quasiparticles near other 



QPTs.l^'26|27j g^jj these cases, the ordered phase pre- 
serves the Ca_v symmetry, and the velocity anisotropy is 
irrelevant so the fixed point is isotropic with emergent 
Lorenz invariance. In the case of a nodal nematic critical 
point, the dispersion anisotropy is enhanced at criticality. 
This has unforeseen effects on the single-particle spectral 
function. 

We conclude with a few observations on the possi- 
ble relevance of a nodal nematic quantum criticality to 
the cuprate physics, recognizing the danger of extrap- 
olating large N results down to = 3^3 The previ- 
ously cited evidence of a nematic phase YBa2Cu30j, near 
y ^ ^.h fronP' transport anisotropy measurements and 
from inelastic neutron scattering in untwined sampled, 
involves studies of materials with unusual purity and 
crystalline perfection. Nematic order couples linearly to 
spatial inhomogeneity and disorder, which thus generi- 
cally causes the ordered phase to be replaced by a glassy 
(domain) phase. A glassy phase with anisotropic do- 
mains has been seen in STM studies of underdoped 
Bi2Sr2CaCu208+i^. It is tempting to identify these 
structures, and the "Fermi arcs" seen in ARPES experi- 
ments on the same materiaP^, with the nematic glass. 
However, this interpretation requires an extrapolation 
from the phase with strong SC order where the theory 
is valid to the non-superconducting pseudogap regime 
where the experiment is carrie d out. M oreover, this in- 
terpretation is far from unique.l^^IIinil A critical test of 
this interpretation is that as the disorder is made increas- 
ingly weak, there will be a crossover from a short range 
disorder behavior with a nodal point that is broadened 
in all direction, to a nematic glass regime with a longer 
T = arc length but a narrower width (with Eq.(5.3l 
as the limiting case). The existence of a nodal nematic 
phas^^, related to ordered or fluctuating stripes, is mod- 
erately cleai!^ in La2-a:Sra;Cu04 and related materials 
near x = 1/8, but it seems likely that the critical phe- 
nomena in this system will be considerably affected by 
quenched disorder. 

The structure of the single particle spectral function 
at the nodal nematic QCP shown in Figj3] with sharp 
qp's within a narrow wedge around the FS but otherwise 
broad, may provide a consistent explanation for both 
the observation of broad qp' s in A RPES— and of sharp 
qp interference peaks in STlvPSEUon Bi2Sr2CaCu208+5. 
Since the wedge is so narrow due to the extreme dis- 
persion anisotropy, it is unlikely to be observed in a di- 
rect momentum space probe. Hence in ARPES, one is 
likely to only observe broad qp's. However the inter- 
pretation of peaks in the Fourier transform of STM on 
Bi2Sr2CaCu208+5 in terms of qp interference invoke the 
notion of coherent qp's at the tips of the equal energy el- 
lipses (or "b anana" s as it is frequently referred to) along 
the FS^SEU. The narrow wedge provides a natural and 
possibly unique mechanism for this interpretation. This 
aspect will be further analyzed in a future publicatiorP^. 
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